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Abstract: We study BPS states in a marginal deformation of super Yang-Mills 
on R x S 3 using a quantum mechanical system of g-commuting matrices. We focus 
mainly on the case where the parameter q is a root of unity, so that the AdS dual 
of the field theory can be associated to an orbifold of AdS 5 x S 5 . We show that in 
the large iV limit, BPS states are described by density distributions of eigenvalues 
and we assign to these distributions a geometrical spacetime interpretation. We go 
beyond BPS configurations by turning on perturbative non-g-commuting excitations. 
Considering states in an appropriate BMN limit, we use a saddle point approximation 
to compute the BMN energy to all perturbative orders in the 't Hooft coupling. We 
also examine some BMN like states that correspond to twisted sector string states 
in the orbifold and we show that our geometrical interpretation of the system is 
consistent with the quantum numbers of the corresponding states under the quantum 
symmetry of the orbifold. 
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1. Introduction 

The AdS/CFT correspondence, in the most celebrated and understood example, 
states the equivalence between the dynamics of type IIB string theory on AdS$ x 
S 5 and M = 4 super Yang-Mills theory in four dimensions Jl|]. According to the 
correspondence, classical gravitational physics should be somehow present in the 
strong 't Hooft coupling regime of the dual gauge field theory. This means that one 
should be able to see the emergence of the AdS§ x S 5 geometry and locality in the 
strong 't Hooft coupling regime and large N limit of Af = 4 super Yang-Mills. 

A first step towards the emergence of locality in the S 5 was taken in |J, by 
studying the dynamics of 1/2 BPS configurations of M = 4 SYM. This problem in 
the field theory side is reduced to a matrix model quantum mechanics for a single 
normal matrix and with first order dynamics. This system is characterized by the 
matrix eigenvalues of this holomorphic matrix, that turns out to be normal, and one is 
led to the study of N free two-dimensional fermions in a lowest Landau level system, 
perturbed by a harmonic oscillator potential. Different 1/2 BPS configurations are 
accounted for by different incompressible fermion droplets in a two-dimensional phase 
space. The ground state of this system, i.e. the circular droplet, is associated to the 
AdSs x S 5 background. In this picture, eigenvalues separated from the droplet and 
holes in the droplet are supposed to represent D-branes in AdS 5 x S 5 . Separating 
several eigenvalues together from the droplet corresponds to placing several D-branes 
on top of each other. The stack of a large number of D-branes should be represented 
by a supergravity solution. 

An impressive and surprising confirmation of this picture was presented in [[|, 
where all regular 1/2 BPS solutions of type IIB supergravity are constructed from 
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some boundary conditions in a two-dimensional section that coincides with the pic- 
ture of fermion droplets in a two-dimensional space. We mentioned that one can 
start to see the emergence of locality in the S 5 in this picture. In this regard, the 
edge of the circular droplet is identified with an equatorial circle of the S 5 , and for 
more general solutions, the plane that describes the fermion droplet configurations 
is a submanifold of the supergravity solution. 

Similar ideas were put forward in but in this case 1/4 and 1/8 BPS con- 
figurations of M = 4 SYM were considered instead. For 1/8 BPS configurations, 
the field theory problem is reduced to a matrix model quantum mechanics for three 
commuting normal matrices. This model is arrived at by realizing that the BPS 
configurations of the field theory on S 3 are related to the moduli space of vacua 
of the field theory. For Af = 4 SYM this means that one has to consider systems 
that satisfy the F-term and D-term constraints, and this leads us to the study of 
six hermitian commuting matrices, that are paired into three normal matrices once 
a choice of a particular M = 1 supersymmetry has been made. 

Again the dynamics is described in terms of the eigenvalues of matrices, but they 
cannot be accounted for as free fermions anymore. Instead one ends up describing 
the system by holomorphic wave functions of N six- dimensional bosons in a har- 
monic oscillator potential, subject to a repulsive interaction. This effective repulsion 
arises from measure terms in the transformation from general commuting matrices 
to diagonal matrices and it is the usual repulsion of eigenvalues that arises from the 
Vandermonde determinant in matrix models. In the large N limit, one can think of 
the system as a thermodynamic system and obtain a density distribution of eigen- 
values in a six-dimensional phase space. For the ground state, the eigenvalues are 
uniformly distributed in a round hollow S 5 . This was argued to correspond exactly 
to the S 5 near horizon factor in the dual geometry. 

If non-BPS modes of finite energy are turned on, they are associated to off- 
diagonal modes of the matrices. These off-diagonal modes to a first approximation 
can be treated perturbatively and they become bi-local in this S 5 geometry. This 
is because they can be drawn as an arrow connecting two different eigenvalues, that 
become associated geometrically with points on the sphere (the sphere depicts the 
distribution of eigenvalues after all). In this setup they can be interpreted as bits 
of massive strings and the S 5 can be directly identified with the S 5 of the geometry 
AdS^ x S 5 that sits at the center of AdS space in global coordinates. This picture can 
be used to successfully reproduce geometrical calculations for some moving strings. 
For instance, the energy of BMN string states [|[] can be obtained in this picture 
by a simple saddle point approximation ||. This formalism for doing calculations 
presents a new viewpoint of how the AdS spacetime geometry is encoded in the dual 
CFT for the particular case of N = 4 SYM on S 3 x R. 

It is natural to ask if this description of the emergence of geometry and gravity 
from SYM is unique to M = 4 SYM or not. In this paper we show that the emergence 
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of gravity from the matrix model describing the BPS configurations of the field theory 
can be extended to other realizations of the AdS/CFT correspondence. What changes 
with respect to M = 4 is that the set of matrices that we need to consider do not 
commute any longer. However, they do not differ significantly from commuting 
matrices, so that one can still talk about distributions of eigenvalues that are given 
by block diagonal sets of matrices describing the scalar fields of the theory. This is 
related to the notions of geometry arising from describing moduli spaces of D-branes 
at Calabi-Yau singularities, as explained in [0. 

Clearly, the simplest examples of such geometries are orbifolds of Af = 4 SYM 
||, and one would want to study the dynamics of branes at these singularities with 
the techniques developed in j3J. 

Some of these orbifolds give very simple theories that have the same field content 
as M = 4 SYM and are natural candidates to investigate. These turn out to be field 
theories on the moduli space of marginal deformations of M = 4 SYM, preserving 
M = 1 supersymmetry ||. These particular deformations of M = 4 SYM have a 
surprisingly rich dynamics and have been extensively studied in the past. The gravity 
duals for those deformations with U(l) 3 global symmetry were recently constructed 
1(J, and this has sparked new interest in the study of these field theories. Also, 



when the deformation parameter q takes special values, the gravity dual is simply 
the near horizon of D-branes in orbifolds with discrete torsion JTT , [12], [13] , which is 
simply AdS$ x S 5 /T, where T is a particular discrete symmetry group Z n x Z n for q 
given as an n-th root of unity. 

We will consider these deformations in the particular case where q is a root of 
unity and we will show the emergence of the orbifolded geometry from a matrix 
model, in a similar vein to the derivation of the geometry of the five sphere in the 
AdS^ x S 5 case. We also obtain within this picture the all loop anomalous dimension 
of the BMN limit of states in the Lunin-Maldacena background and the orbifold 
geometries, plus we show that some peculiarities of the saddle point evaluation of 
energies of these states are required for consistency of the geometric interpretation 
as string states in AdS$ x S 5 /T with some fixed quantum numbers with respect to 
the quantum symmetry of the orbifold. 

The paper is organized as follows. In section we derive a matrix model of 
g-commuting matrices describing the BPS configurations of a special deformation of 
M = 4 SYM and show the emergence of the orbifolded S 5 as a density distribution of 
eigenvalues. In section ^| we show that non-BPS modes in this setup are localized in 
the orbifolded S 5 and interpreted as massive string bits. In section |] we compute the 
energy of BMN states to all perturbative orders, and find that similar states that are 
in the twisted sector correspond to certain non-trivial loops in the orbifold geometry 
that lift to open paths in the covering space. In section [5] we discuss and summarize 
the results of this paper. 
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2. Quantum mechanics for g-commuting matrices 

To begin, we will study a one parameter set of marginal deformations of M = 4 SYM 
that has U(l) 3 global symmetry. These were originally shown to be conformal by 
Leigh and Strassler || . These depend on a complex parameter q, that we will further 
specialize to be a root of unity, and they present the following superpotential for the 
three adjoint superfields X, Y and Z, 

W q = tr(XYZ) - qtr(XZY) . (2.1) 

When the deformation parameter q = r exp(2irif]) takes the special values r = 1 
and f3 rational, the CFT is related to geometries of orbifolds with discrete torsion 
[ |Tl"| , m , [TB| ] . For simplicity we will focus exactly on these deformations such that q is a 
primitive n-th root of unity, and we will try to see how our calculations depend on q. 
Our goal is to be able to say that the ground state of the CFT dual of these theories is 
given geometrically by AdS^, x S 5 /T and that we can reproduce part of the spectrum 
of strings on this geometry from our considerations. Our ultimate goal is to show 
that the techniques developed in £|, § can be applied in a broad class of examples 
and are not restricted to having a quantum system with 32 supersymmetries. 

The line of thought explored in [|J begins by considering the chiral ring of the field 
theory on Euclidean space, and trying to understand what states are associated to 
the operators in the chiral ring under the operator state correspondence. Chiral ring 
operators preserve some of the supersymmetries, so the dual states are BPS states. 
Therefore, in the AdS/CFT correspondence we want to consider BPS configurations 
of the field theory on S 3 , that are dual to the set of chiral primary operators. These 
operators, in the free field limit can only be built from the constant modes of the 
fields on the S 3 . One can additionally make some other such states by turning on 
some of the fermion fields. However these are zero in a naive semiclassical treatment, 
so we will ignore them for the purposes of this paper. From these considerations we 
get an effective dimensional reduction to a matrix model whose degrees of freedom 
are the constant scalar modes of the fields on the S 3 . Focusing on these modes alone, 
we get the following effective action 

S = ti 2 J dttr (\D t X\ 2 + \D t Y\ 2 + \D t Z\ 2 - \X\ 2 - \Y\ 2 - \Z\ 2 - V D - V F ) , (2.2) 

where Vd and V F are the potentials resulting from integrating out the auxiliary D- 
terms and F-terms in the M = 1 lagrangian. The factors of 7r result from integration 
over the volume of the S 3 . Moreover, the BPS constraint is of the form 

A = J (2.3) 

where A is the Hamiltonian of the field theory (and also the generator of dilatations 
of the Euclidean field theory under the operator-state correspondence map), while J 
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is the generator of the R charge with a slightly different normalization that matches 
the BMN conventions ||. In the free field limit we satisfy this constraint by taking 
positive frequency solutions of the holomorphic fields of the M = 1 theory and in 
this limit we ignore Vp, and Vp terms in the action. With the normalizations above 
we get that the solutions have the following time dependence 

X(t) = X(0)exp(it) (2.4) 
Y(t) = Y(0) exp(tt) (2.5) 
Z(t) = Z(0) exp(it) (2.6) 



To consider classical solutions of the BPS constraints in the interacting case, we 
notice that the solutions above satisfy A = J if one also imposes the constraints 
Vp = Vp = 0. In the usual J\f = 1 supersymmetric field theories, this corresponds 
to the classical moduli space problem of finding the supersymmetric vacua of the 
supersymmetric field theory. For theories associated to D-brane actions, this can be 
solved in two steps. First solve the F-terms, that are given by matrix equations, and 
afterwards on the set of these solutions one does complexified gauge transformations 
to try solve the D-term constraints. 

For this particular field theory, the characterization of the moduli space of su- 



persymmetric vacua of these marginal deformations have been studied in [13], [14 



The classical moduli space is determined by solutions of the F-terms constraints, up 
to GL(N,C) equivalence (the complexification of the gauge group, as is familiar in 
supersymmetric theories) 



XY - qYX = , 
YZ - qZY = , 
ZX - qXZ = . 



(2.7) 



This means that we have to look for NxN g-commuting matrices. We are considering 
the special case q = exp(27rz/n). For this case, an n-dimensional representation of 
the algebra ( |2.7| ) can be obtained in terms of the clock and shift matrices P and Q 



P 



/10. 
q . 








Q 



\0 ... q n ~ x ) 



/0 ... 1\ 
1 ... 

\0 10/ 



(2.1 



The algebra ( |2.7D is solved by the following representation: simply take X ~ P, 
Y ~ Q and Z ~ P n ~ 1 Q n ~ 1 , and this is an n-dimensional irreducible representation 
of the algebra. Noticing that X n ,Y n , Z n ,XYZ are proportional to the identity 
on these representations, one can identify these as coordinates of a commutative 
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geometry, that turns out to be exactly the ring of holomorphic functions on the 
orbifold C 3 /Z n x Z n . 

There are also branches of fractional branes where two out of the three matrices 



are zero X = Y = lets say |13], |14| . These can have more moduli associated to 
them as the eigenvalues of Z become unrelated to each other 1 . 

However, it is easy to show that their orbits in the configuration space of matrices 
(ignoring the equivalence of representations under GL(N,C)) always have higher 
codimension than the matrices shown above. Since we are trying to understand the 
ground state of the matrix model we discussed, one can argue that these branches are 
less important because they have higher codimension on the set of allowed matrix 
configurations. We will ignore these branches in what follows, assuming that for 
the ground state wave function, the other branch made only of branes in the bulk 
dominates. This will be shown to be self-consistent later on in the paper. Moreover, 
these extra fractional brane branches are only present for singular three dimensional 
CY spaces where the singularity at which the branes are placed is not isolated, so 
they are not relevant for isolated singularities. 

Let us now also assume for simplicity that the rank of the gauge group is N = mn. 
In this case, it is straightforward to construct an ^-dimensional representation of 



the same algebra, with tensor products of matrices |]15 



X= -= A X ®P, (2.9) 



v 7 ^ 



n 



Y = ^==A 2 ®Q, (2.10) 
V In n 



1 



V2^ 



A 3 <g> P n ~ x Q n ~ x , (2.11) 



n 



where A 1 , A 2 and A 3 are three commuting holomorphic m x m matrices. These 
solutions are just direct sums of the irreducible representations found above, and 
they solve the F-term equations automatically. Indeed, one can show that the generic 
representation of m D-branes in the bulk away from the singularities of the geometry 
corresponds exactly with such a representation, and that any other general solution 
can be transformed into this form by a gauge transformation. This is explained in 
detail in [|TB|, p]. Here we will just use these results. In [|TB|, [7| was also noted that 
there are special solutions of representations where two of X, Y, Z are zero and the 
third is a c-number. These correspond to fractional branes at the singular locus. 
These solutions are ignored above. 

The normalization factor above for the parametrization of the matrices was in- 
cluded for later convenience (this normalization guarantees that the action for the 



1 These branches have been studied non-perturbatively in the work of Dorey and collaborators, 
see [l^, [I?]] and references therein. 
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matrices A* has a canonical normalization when the theory is compactified on the 
sphere) . 

Now, the moduli space of supersymmetric vacua is parametrized by 3m eigen- 
values of the commuting matrices A 1 , A 2 and A 3 . We should notice that, if we call this 
representation ^(A 1 , A 2 , A 3 ), representations RiqA 1 , g _1 A 2 , A 3 ) and RiqA 1 , A 2 , g _1 A 3 ) 
are equivalent under similarity transformations. This is, these configurations are 
equivalent under gauge transformations. Imposing this equivalence, the eigenvalues 
parametrizing the moduli space take values in the orbifold C 3 /Z n x Z n . This proce- 
dure can be implemented at the level of wave functions for the matrices by requiring 
that the wave function is invariant under these identifications of configurations. 

We can now obtain a matrix model from the scalar sector of the M = 1 marginal 
deformation of M = 4 SYM action in question. We define the field theory in R x S 3 
and study BPS configurations keeping only the s-waves modes of the g-commuting 
scalar fields ( |2.10|) . This follows the techniques used in |4[], where these classical 



states satisfying the F-term constraints on the S 3 are related to operators on the 
chiral ring of the field theory. 

To obtain the matrix model we write an ansatz with spherical symmetry for the 
fields that respect these F-term constraints into the effective action of the field theory 
on S 3 . For solutions with this ansatz, that moreover satisfy the D-term constraints as 
well, the D-term and F-term potential term in the SYM action on S 3 vanishes. The 
solutions as written above satisfy all the necessary constraints, provided that the A* 
matrices are normal matrices and that they commute with each other. These fields 
are still massive in the reduced matrix model because these scalars are conformally 
coupled to the round metric on the S 3 . We then perform the integral over the S 3 
and we get a quantum mechanics matrix model with a quadratic potential, 

S = 2vr 2 J dt (\D t X\ 2 + \D t Y\ 2 + \D t Z\ 2 - \X\ 2 - \Y\ 2 - \Z\ 2 ) 

= J dt itr (| AA 1 ! 2 + | AA 2 | 2 + | AA 3 | 2 - (A 1 ) 2 - |A 2 | 2 - |A 3 | 2 ) . (2.12) 

The normalization of the matrices cancels factors of n from the trace of the identity on 
the set of P, Q matrices. Since the A* matrices are mutually commuting, and normal, 
we can diagonalize them simultaneously with a unitary transformation, together with 

A\ 

This is a gauge transformation on the set of configurations. This reduces the 
degrees of freedom from generic commuting matrices to their eigenvalues. Classically, 
the system is reduced to m decoupled harmonic oscillators in three dimensions. 

However, the emergence of gravity as detailed in j| takes place in the quantum 
mechanical matrix model. The reduction of degrees of freedom to the eigenvalues 
of the A* matrices is valid only if we include the measure factors coming from the 
volume of the gauge orbit. Because of this measure, the harmonic oscillators will 
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no longer be decoupled. Since there are global transformations that permute the 
eigenvalues, the wave functions also have to be symmetric under these permutations. 
As done in |3J] , this system will be interpreted as a set of m interacting bosons in a 6 
dimensional phase space. This requires looking at the BPS constraint, that only keeps 
the positive frequency modes of the A* matrices. This procedure effectively converts 
the classical moduli space into a symplectic geometry, because for these solutions 
^xy.z ~ A, F, Z and one has an effective magnetic field on the moduli space that 
makes the problem similar to a quantum hall system: we want to quantize only the 
lowest Landau level problem. 

Describing the g-commuting matrices as fl2.10|) with diagonal A* is a gauge choice. 



The quantum modifications of the free Hamiltonian that we need arise from comput- 
ing measure effects from changes of variables. In the original variables, the Laplacian 
of the Hamiltonian has a cartesian coordinate form. Taking into consideration mea- 
sure effects is very similar to the problem of writing the Laplacian of Euclidean flat 
space in spherical coordinates, using the integration measure of spherical coordinates. 

The measure terms we are after come from calculating the volume of the gauge 
orbits in question. To compute this volume we set matrices in the form ( p.lC| ) 
with diagonal A 1 and perform an infinitesimal gauge transformation with the broken 
generators, i.e. transformations that turn on off-diagonal components in matrices A 
and give rise to general combinations P a Q 13 in the tensor products 2 . These variations 
are 



■ 



5x = q - p e a J(xl - /A*)r a6 ® p« + V 

5Z = gP&S(Xl - q a ^Xl)T ab ® P«- V 1 , (2.13) 

where O^f parametrize the transformations generated by T ab ® P a Q lS and T ab is a 
m x m matrix whose components are (T a b) cl i = 5 ac 5bd- 

As we said, the vectors of complex eigenvalues A a = (A*, X 2 , A^) take values in 
the complex orbifold C 3 /Z n x Z n because of gauge identifications between different 
configurations. Roughly, we could plot all these eigenvalues in a n 2 -th wedge of the 
whole C 3 . For expressing the measure, it is convenient to distinguish among the n 2 
images in the complete C 3 . So, we define A"" 5 = (A^g a+/3 , A^g" a , A 3 g _/3 ). The volume 
associated to the variation is equal to the length of the vector (A a — \b) a ' 13 '■ Then, 
the measure turns out to be 

^^n'n^-^i-nn^-^T- (2.14) 

o,/3 a a, /3 a<b 

The prime in the first product is to remark that the case a = (3 = n is excluded. 

At this point, we can define new indices A,B,... running from 1 to mn 2 . We also 

2 To refer to different kinds of indices we use different sets of letters. Indices a, b, . . . = 1, . . . , m 
and indices a,0,... — 1, . . . , n. 
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call Xa the mn 2 vectors A"'". Moreover, the measure (|2. 14 ) resembles the vectorial 



generalization of the Vandermonde determinant that appeared in |4]], evaluated for 
the mn 2 vectors A^. However, many |A^ — A#| 2 are missing. In the complete Van- 
dermonde determinant factors are repeated because only m of the mn 2 vectors A"'^ 
are arbitrary, while the rest are the images under the action of Z n x Z n . Measure 
( |2.14D does not present such a repetition of factors and coincides exactly with the 



n 2 -th root of this determinant, where each different length is counted only once 

f= [\{^a-^bA ■ (2.15) 

\A<B J 

We could have guessed this form of the measure by using the method of images on 
the covering space and taking the n 2 -th root of the result to avoid over-counting of 
states. 

Now, we use this measure factor to define the reduced Hamiltonian in the eigen- 
value basis 

H = J2-—2 V ^ a + ^\\ a \ 2 . (2.16) 

a 

We will look for the ground state wave function of Hamiltonian ( |2.16|) and give a 
probabilistic interpretation of its square modulus as a Boltzman distribution in the 
limit of large m. As in the undeformed case || , we can show that the ground state 
for m decoupled harmonic oscillators, 

^ ~exp(-^|A a | 2 /2), (2.17) 



is an exact eigenfunction of (|2.16|) that has the correct symmetry under exchange of 
particles and also under the discrete group orbit V . In what follow, and to simplify 
the notation, we write the vectors of 3 complex dimensions A a as vectors of 6 real 
dimensions x a . Moreover, the vector gradient V a now refers to the gradient associated 
to the 6 real coordinates. Using that V a ipo = —x a 4>o, we have 

a ~ 

= Yl ^(£aV a/ u 2 )Vo + 6m^ • (2.18) 

As in the undeformed case, following 0], the measure is a homogeneous function, 
now of degree m(mn 2 — 1). Then, /i 2 is an eigenfunction of the operator x a V a and 
i/jq is an eigenfunction of the reduced Hamiltonian. Since, is rea l an d positive one 
expects this eigenfunction to be the ground state of the system. The states dual to 
primary fields should arise from multiplying the above ground state by holomorphic 
polynomials in the X, Y, Z fields, subject to the moduli space constraints. This 
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effectively results in a subset of the polynomials of the A*, because the traces also 
include traces over the P, Q matrices part. These vanish unless the corresponding 
trace is proportional to the identity. 

Orthogonality relations with another eigenfunction and expectation values of 
observables have to be computed using the integration measure p? we have found. 
For that reason, it is convenient to absorb a p factor into the definition of the wave 
function, 

$ = fiip . (2.19) 

Then, after this similarity transformation, we get the usual integration measure 
Q(i 6 x a associated to wave functions ip. The factor p is also symmetric under the 
exchange of all the vectors x a and under the discrete identifications of eigenvalues. 
So the particles are also m identical bosons with respect to the new measure. To 
consider the emergence of geometry, we are instructed to place the m bosons in the 
same phase space and study the density distributions of particles for the given wave- 
functions described by the ip. We can consider the modulus square of a given wave 
function as a probability distribution of the m bosons in the 6 dimensional space 
C 3 /Z n x Z n . For the wave function ip we have 

-a,P\ 



x b 



|4 2 | ~ exp - l^| 2 - J2 lo S l f « - ^\ + 2 J2 E lo S l f « - 

V a \ oi,P J a<b a, (3 

(2.20) 

The computation of position observables from the square of the wave function is 
equivalent to a statistical mechanics problem of calculating the partition function 
for m bosons in an external quadratic potential. These bosons are subject to a 
logarithmic repulsion between them and with respect to their images through the 
action of Z„ x Z„. We consider the thermodynamic limit m — > oo, where the bosons 
form some continuous and positive distribution density p on the phase space of a 
single particle and look for the most probable p. In this continuous limit, sums over 
a are converted into integrals over C 3 /Z n x Z n , or integrals over C 3 if we take care 
of normalizing the density at a point and its images properly. 

£(...) - / d e xp(x)(...) (2.21) 

a Jc 3 /Z n xZ n 

m = d 6 x p{x) (2.22) 

ic 3 /Z n xZ n 

The probability distribution in terms of the density p is 

\ipl\ ~exp(- / d 6 xp(x)x 2 + [ d & xd 6 yp{x)p{y)^\og\x-y a ^\ J . (2.23) 

Maximizing the probability density ( |2.23| ) we can find the most likely boson distri- 
bution, exactly as in H. Varying the argument of the exponential we obtain the 
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following integral equation for p(x) 

x 2 + C = 2 f d 6 y po(?/)^log \x a - Ta P \ • (2.24) 

The constant C is a Lagrange multiplier enforcing the constraint in the total number 



of bosons (|2.22j ). In 6 dimensions the function log |af — y\ is proportional to the 
Green's function for the operator (V 2 ) 3 . If we operate naively on both sides of Q2.24| ) 
with this operator we obtain p = . From this expression we conclude that the 
distribution po has singular support, because the naive manipulations are only valid 
in the case that p is differentiable. In the covering space of the orbifold, we see 
a spherically symmetric Boltzman gas of particles, identical in nature to that one 
found in j4[]. 

Because of symmetry, the simplest ansatz for po we can make is a singular spher- 
ically symmetric distribution at radius r on C 3 /Z n x Z n . So, the eigenvalues of the 
ground state are uniformly distributed in an orbifolded 5-sphere S 5 /Z n x Z n , 

P0 ~ m rlVo\(SyZ n x Z„) ' (2 ' 25) 



where the normalization was chosen so that the constraint Q2.22p is fulfilled. It is 



expected that this is the unique saddle point for the distribution of the density of 
eigenvalues (this has not been proved yet for this problem). 

This S" 5 /Z n x Z n can be identified with the compact factor of AdS§ x 5' 5 /Z n x Z n , 
which is the geometrical background of the string theory dual to the g-deformation 



of SYM we are considering [13, 14]. In the next section, we will show how the energy 
of simple non-BPS excitations are localized in the S' 5 /Z ra x Z ra and we will be able 
to interpret them in terms of bits of massive strings. 

Now we substitute this ansatz into the argument of the exponential (we call it 
7i) in ( 2.23| ), and minimize with respect to r . Since the distribution is supported 



at ro and uniformly distributed in the angular variables, it is easy to isolate the r$ 
dependence. We obtain, 

H = mr 2 — m 2 n 2 logr + k , (2.26) 

where k is a very complicated integral on the angles of the S" 5 /Z n x Z n but inde- 
pendent of r . Now, minimizing 7i with respect to r we easily find the critical 
radius 

mn 2 nN , n 

r " = V — = v — ■ (2 - 27) 

This result is relevant for the computation of BMN excitation energies we will carry 
through in section ^. At first sight it seems that this result depends on n. However, 
this is an artifact of the normalization factor of n used in ( p. 11 ). If we restore the 
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standard matrix normalization, we find that the radius of this sphere is independent 
of n (this translates into the radius being independent of q, for q unitary). However, 
some of the intermediate steps to derive the result would involve factors of n in various 
other places (the Hamiltonian and the wave function for example). This result is 
crucial for us. This is what ensures that calculations can be analytic functions of 
q. For analytic functions, knowing the rational values of the function lets us find 
the result for arbitrary values of the phase of q by continuity. This is what will let 
us extend some of our analysis to the non-BPS case and to explore states in the 
Lunin-Maldacena background [[lO|j . 

Most of the other BPS wave functions of the system are given by holomorphic 
functions of the A multiplying ipo that satisfy the boson statistics and that are in- 
variant under the identifications produced by T. These symmetries are all required 
from the residual gauge identification of configurations. In general, one obtains in 
this way the chiral ring associated to gravity modes in the bulk of the spacetime 
(untwisted states in the chiral ring). These can also be interpreted as wave functions 
in the covering space. If one follows the gravity interpretation of these states in the 
AdSz, x S 5 case, one finds that coherent states of single trace operators correspond to 
deformations of the S 5 geometry that respect the discrete symmetries of the orbifold, 
and hence they correspond to supergravity solutions that deform the shape of the 
S 5 /T. 

To include also the twisted sector BPS states, one also needs to consider having 
some small fraction of the eigenvalues to be located in the fractional brane branches. 
Understanding exactly how this works is beyond the scope of the present paper, 
but it is a very interesting problem to consider, as these deformations smooth out 
the singularities of the S 5 /T space. Some progress for solutions with branes on 
the Coulomb branch of the field theory has been made at the level of supergravity 



solutions with Lunin-Maldacena boundary conditions in [191, but these preserve 



many symmetries that make the configurations non-generic. 
3. String bits 

In this section, we want to describe the emergence of massive stringy modes in the 



CFT as in ||, see also pC§ for an alternative viewpoint on some of these computa- 
tions. The first thing to do is to go beyond BPS configurations. So, let us include in 
the action the potential D-terms and F-terms and consider situations where we are 
in a near BPS state. The action for the s-wave scalar modes is 



S = n 2 



J dtti (\D t X\ 2 + \D t Y\ 2 + \D t Z\ 2 - \X\ 2 - \Y\ 2 - \Z\ 2 - V D - V F ) ,(3.1) 



-12- 



where 



V D = ^tr [\[X,X] + [Y,Y] + [Z,Z}\ 2 ) , (3.2) 
V F = ^i X (\[X,Y] 9 \* + \\Y,Z\ 9 \> + \[Z,X\ q \ 2 ) . (3.3) 

Being near BPS, means that we can attempt to treat the degrees of freedom that take 
us away from being BPS as a perturbation. These should be given by infinitesimal 
deformations of BPS configurations. These configurations result from turning on 
the other matrix components of the fields X, Y, Z that can not be turned on in a 
BPS configuration. In perturbation theory we treat these non-BPS modes as fast 
degrees of freedom attached to pairs of eigenvalues. For a non-BPS component in 
the matrix position a, b, one can associate the infinitesimal matrix component to the 
eigenvalues of the matrix in the diagonal components a and b. The eigenvalues in 
this setup are treated as slow degrees of freedom. As a first approximation, we treat 
the off-diagonal modes as free fields coupled to a background of eigenvalues, and we 
ignore the back-reaction of the eigenvalues to the string bits. 

In order turn on a non-BPS mode we should include in at least one of the holo- 
morphic fields X,Y,Z an off-diagonal component in the matrices A 1 , A 2 , A 3 and/or 
a general P a Q@ matrix in the second factor of the tensor product. Let us consider, 
for instance, 

X' = X + 5X$g (3.4) 
where X is the same as in ( |2.1(J| ) and the non-BPS piece is 

sx$$ = -i= ® P 1+a Q" , (3-5) 

V 2n z n 

where the matrix A^ Q/3 has components (A^ b at3 ) c d = Xl' b af3 5 a c5bd- This mode is off- 
diagonal, unless a = b and a = f3 = n. It is considered to be off-diagonal with 
respect to the eigenvalue decomposition of X if the P, Q matrix component does not 
commute with X, or does not g-commute with Y and Z according to the F-terms 
of the field theory. Because in general the modes 5X, 5Y, 5Z can mix, the problem 
of finding the correct modes can be complicated. However, we can use the fact that 
this theory is obtained from an orbifold of the M = 4 SYM theory to diagonalize 
the problem. It turns out that after mixing is taken into account, the Hamiltonian 
associated to these non-BPS modes is 

^r=^inri 2 +^(-r) 2 iAri 2 , (3.6) 

where the frequency of the mode turns out to be 



(3.7) 



In the end, we have to evaluate these energies using the distribution of eigenvalues 
obtained in the previous section. The n 2 non-BPS modes for two given points in the 
orbifold can be seen as the n 2 possible straight lines connecting the first point with 
the n 2 images of the second. This can be easily visualized in the cone C 2 /Z 2 . In the 
Figure [I] we draw the two possible straight lines connecting two points on the cone. 





Figure 1: Straight lines between two points on C 2 /Z,2- 

The interpretation of these modes as bits of massive strings goes exactly as in 
the undeformed case M. According to the AdS/CFT correspondence, in the limit 
of large 't Hooft coupling A, the system should be described in terms of classical 
gravity with some string probes on the background. In this regime, using normalized 
distances r, the energies of the non-BPS modes are order 



E 



a/3 



A 1/2 lr„ - 



-o.,a—j3 1 



(3.8) 



where we have used the radius normalization found in the previous section to rewrite 
the mass terms in terms of the 't Hooft coupling, and unit vectors on the sphere 
r a . At strong coupling, the distance between the points dominates the energy of 
the string bits. Finite small energy string bits (energy of order one) require very 
small \r a _ r ~ a,a_ ^| anc l the segments between r a and r^" ' 0- ^ become tangent to the 
5-sphere. These segments have to be joined into closed loops because we need to 
make a gauge invariant state. 

We can consider these small segments as short bits of strings localized on the 
geometrical S 5 /Z* n x Z n . In the case of strings on AdS^ x S 5 , it was shown in |§ that 
this description in terms of free short string bits was useful for the BMN limit, but 
longer string bits should be treated as interacting fields. Even though long string 
bits are interacting, the leading order perturbative evaluation of their energies gives 
the correct order of magnitude estimates for string energies and one can also use 
them to get some intuition on the eigenvalue geometry they are associated with. 



4. Energy for BMN states 

In this section we repeat a similar analysis to that of H] to obtain an all-loop expres- 
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sion for the BMN operators energy in the g-deformed theory. To do that, we compute 
the approximate energy of the state resulting when we excite the ground state Q2.17Q 
with a BMN or BMN-like operator. BMN operators in marginal deformations of 
M = 4 SYM have been already studied in the literature [p3L p4L 25 1, by considering 



appropriate deformations of the usual spin chain describing the dilatation operator. 
For instance, consider the dilatation operator in the two spin sector, at the 1-loop 
approximation, and with a deformation parameter q = rexp(2iri(3). For a real de- 
formation q = r, the dilatation operator turns out to be a periodic XXZ spin chain 
[f23| , |24| . The anisotropy parameter is a function of the deformation parameter such 
that A(r) < 1 and then the spin chain is always ferromagnetic. When the defor- 
mation parameter presents a non-trivial phase, it is possible to perform a position 
dependent change of basis. In this particular new basis the Hamiltonian of an XXZ 
spin chain is recovered and the phase of the deformation is totally encoded in some 
twisted boundary conditions and a deformation in the cyclicity condition pi ]. More- 



over, in the case we are interested in, i.e. when the deformation parameter is a pure 
phase q = exp(27ri/3), the chain becomes an XXX spin chain with twisted boundary 
conditions. The twisting factor is q J where J is the length of the spin chain. 

It turns out that BMN-like operators with two different holomorphic impurities 
are essentially identical to those of the undeformed theory and the momentum 
numbers in the 1-loop energy are shifted kj J — ► kj J ± (3. A consistent BMN limit 
for this spin chain requires that k± f3J should be kept finite in the J — > oo limit. In 
our case, it means that k ± J/n should be finite. So, we will consider the following 
BMN-like operator for exciting the ground state ( p. 17] ), 

j 

O k ~ J2 e M i( fil)^(Z l X Z J ~ l Y) , (4.1) 
z=o 

where (p is the usual BMN phase (p = We take Z to be of the form (|2.11|) and the 
corresponding eigenvalues given by the distribution obtained in section 3. On the 
other hand, fields Y and X include non-BPS modes, that can be treated as creation 

2 a 
v a6 



operators in a quantum system like (|3.6| ). We will refer to them as A^ b a/3 ^ and A^. a ' 



They satisfy the usual creation-annihilation algebra 

[>g,X$] = 80*5^6*. (4.2) 
The above operator is represented by the following state 

(=0 a,/3 a,b 

(4.3) 

The integer s is the reduction modulo n of the number of sites J. We have written 
the wave function tp in the coordinate basis and the non-BPS modes as creation 
operators acting on the vacuum |0) of non-BPS modes. 
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We compute the energy of the above state, as an expectation value, 

E ~ — l — (4 4 

The Hamiltonian (|3.6| ) tells us that the oscillators carry an energy 

rna,l3 l.a/3 . 2.r-a-l.r-/3-l / , r\ 

E ab = UJ ab +^ba ■ ( 4 '5) 

Including the energy of the diagonal piece the total energy is 

E total = J+ (goscj _ ^ 

Then, we have to compute the average energy of the oscillators ( |4.5| ) for the wave 
function ( [4.3| ). This results in the following integral 

/ rposc\ _ ' ' ' M V J_ 

jn^ c |Vio| 2 E aA ^iE^^^ +1 )(A3) ; (A3)^|2 

(4.7) 

As it was done in the undeformed case || , we can compute this average by a saddle 
point approximation. The first thing to notice is that both integral, in the denom- 
inator and in the numerator are dominated by the configurations that maximize 
|-?/>o| 2 - This means that A a and A& should be located exactly on the orbifolded sphere 
= S' 5 /Z n xZ n . More precisely, in the thermodynamic limit sums E a are converted 
into integrals over the complex orbifold C 3 /Z„, x Z n which, using the distribution 
Po{x), are reduced to integrals over fl 5 . 



(E osc ) ~ 



Using A 3 = r cos6'exp(i0), the square of the sum can be rephrased as 



^ e tf V ^(a-/9-l)( A 3y (A 3y-i 



^5Z(cos^)' +/ '(cos^) 2J ^' x (4.9) 



.2./ 

I I 



exp 



z(Z - ( p + - 0' + —(a - p + 1) 



The BMN limit requires J —>■ oo. This serves for improving the saddle point 
approximation thanks to the extra powers of cos 9 and cos 9' in the angular integrals. 
They are maximized when |A 3 | and |A 3 '| take their maximum value on f2 5 , i.e. when 
cos 9 = cos 9' = 1. This corresponds to a geometrical localization for the string with 
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large angular momentum J around a null geodesic. Because of this localization the 
distances in the frequencies of the non-BPS modes are simplified. For instance, 



\x-x> a,a P \ 2 -> |A 3 -/- Q A 3 '| 2 = 4r 2 sin 2 Q (0 - 0' + 2vr(a - . (4.10) 

A similar simplification takes place for the other frequency, leading to the same 
value ( f4.10|) . Moreover, in the large J limit, we can approximate the sum over 
relative phases in (|4.9| ) by a delta function and compute easily the remaining angular 
integrals. The result is that the amount <ft — <$' + 2Tc(a — f3)/n is sharply peaked at 
the value <p + 2n/n, and the energy of the oscillators is 



Using the value we obtained r = \J in the section 3., ip = ^ and taking the 
large J limit, the energy for the oscillators reduces to 



k V 2 



(E osc ) = 2\ll + g 2 YMN{j + -) • (4.12) 

which is exactly the energy of BMN string excitations, with the appropriate nor- 
malization g\ M = 4:7rg s . From the gauge theory side it corresponds to an all orders 
result in perturbation theory. This all loop gauge theory result has been already 
reported in To obtain this result the authors of fl2l| , extended the arguments 



used in p2| for the undeformed case. This includes using the equations of motion, 
but it is not clear if contact terms could spoil the use of these relations. Also, string 
states in the Lunin-Maldacena background have been considered by many authors 
2|, H, |7|, |8|, |9| 



We can also study in the same approximation BMN like states that are not 
reduced to the BMN limit. For example, we can consider the operator tr (Z J Y). 



These operators have also been considered in |3T| (our results agree with those 
calculations). It is easy to see that this trace is zero on BPS states. In this case, 
we interpret Y as the off-diagonal mode. Doing the same analysis as above, we also 
localize on the circle where \z\ = 1, and it turns out that the angle between the 
two eigenvalues that Y connects is exactly given by the argument of q. This means 
that in the covering space, the mode Y goes between an eigenvalue of Az and a 
particular one of it's images. At first sight one might think that one is violating 
the Gauss' law for the corresponding state, because it seems as if the string state 
is an open string. However, this is a non-trivial loop in the orbifold geometry, so 
the correct interpretation is that the associated closed string state is actually in 
the twisted sector of the orbifold. This is exactly as should be. Twisted states of 
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an abelian orbifold have non-trivial quantum symmetry quantum numbers. These 
discrete quantum numbers of the states can be calculated counting powers of X, Y, Z 



modulo n [13|, |I4| . It is easy to see that these jump assignments between images 
are all consistent with the quantum symmetry quantum numbers of states in these 
papers. Moreover, we can ask now under what conditions do some of these states 
survive in the BMN limit. It's easy to show that if we take the BMN limit along Z, we 
need the state to be built only of small bits, and that forces us to have j^Y — j^X = 
mod (n). This reproduces the BMN set of states of the orbifold C 2 /Z n , exactly as 
one would have expected from spin chain model considerations [24|. 



5. Discussion 

In this paper we studied BPS configurations in some marginal deformations of SYM 
preserving M = 1 supersymmetry. We focused our attention to the case in which the 
commutator in the superpotential is deformed to a g-commutator, with q being an n- 
th root of unity. In those cases, the gauge theory is dual to the near horizon geometry 
of N D-branes in the complex orbifold C 3 /Z„ x Z n . The moduli space is described 
in terms of quantum mechanical g-commuting matrices. In the large N limit, this 
is reduced to the study of density distributions of particles (interacting bosons) in 
a 6-dimensional phase space, that results from studying the moduli space of vacua 
of D-branes in the bulk. The ground state corresponds to a spherical distribution of 
bosons in C 3 /Z n x Z n and we identified this distribution as the compact factor of 
the dual geometry. Moreover, we considered non-BPS configurations by turning on 
non-trivial g-commutators. We argued that, as in the undeformed pure M = 4 SYM 
case, the energy of non-BPS modes is localized in this density distribution identified 
with the compact factor of the geometry, and they are considered as bits of massive 
strings. Thus, one of the accomplishment of this paper is that we have an explicit 
example of the extension of qualitative picture jl for the appearance of spacetime 
geometry and locality in the strong 't Hooft coupling regime beyond the standard 
AdS 5 x S r °/Af = 4 SYM. 

Then, we proceeded as in to compute the BMN energies to all orders in 
perturbation theory. The agreement with the string theory computation is complete 
where it is expected. This computation shows that the geometrical interpretation of 
the density distribution of boson is not just a qualitative picture, but that it can be 
used to reproduce concrete string theory calculations. We want to emphasize that the 
approximations we made are good only in the strict BMN limit. If we were interested 
in near-BMN corrections, we would have to deal with string bits of finite length and 
their interactions. The energies of these string states are analytic functions of the 
phase of q (for q unitary), so we can extend them to non-rational values of q and we 
get agreement with other calculation of the corresponding energies. 
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Moreover, we have seen that the geometric picture given by this distribution of 
eigenvalues accounts correctly for discrete quantum numbers of string states, related 
to being in the twisted sector of the orbifold. This is a necessary step to claim 
that we understand the origin of geometry in a string theory setup that is not just 
super gravity. 

Obviously, since we have succeeded in showing that for certain orbifolds of M = 4 
SYM theory we are able to find the correct non-spherical horizon of the geometry, 
it is easy to conjecture that it will work the same way for all such supersymmetric 
orbifolds. Some recent progress in the better understanding of the structure of the 
field theory at orbifolds has been made recently in |3^, [33[ It would be interesting to 
show that this procedure works for other cases as well, such as the field theory ob- 
tained by placing D3-branes at a conifold singularity []34|| . There is some preliminary 
evidence that this is possible in a large class of examples p5| . 
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